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1 Introduction 



The possibility that space-time has more than four dimensions has been considered for many 
years. The Kaluza-Klein approach, widely used in superstring constructions and brane world 
models, relies on the assumption that each four dimensional space-time point is replaced by 
a compact internal manifold. Another possibility which was considered recently is to have an 
infinite fifth dimension with a negative cosmological constant, the four dimensional universe being 
confined to a 3-brane [1] on which gravity is localised. A generic feature of such models with 
continuous extra-dimensions is that the field excitations in the internal manifold are seen by four 
dimensional observers as a tower of infinitely many massive particles. This holds in particular 
for excitations of the metric with non zero momentum in the internal dimensions, resulting in 
a tower of massive gravitons. On the other hand, theories with a single (or a finite number of) 
massive graviton have also been considered. In particular it is well known, that there is only one 
ghost-free quadratic action for a massive spin 2 field, first introduced by Pauli and Fierz [2] . Going 
beyond quadratic order, one is lead to theories with several coupled metrics [3, 4]. Such theories, 
which we will call generically here multigravity, have been investigated in the past in relation to 
strong interaction [3], and more recently by various authors [5, 6] in particular in relation with 
discrete extra-dimensions [7, 6] . Discrete extra dimensions have also been considered in the spirit 
of non-commutative geometry [8] . 

Theories with massive graviton excitations are likely to suffer in general from several patholo- 
gies which could prevent them from being very useful to describe the real world. It includes the 
celebrated van Dam-Veltman-Zakharov (vDVZ) discontinuity [9]^, as well as the propagation of 
ghost-like fields when one considers an action which goes beyond quadratic orders [4]. Nonethe- 
less some of these theories have promising phcnomenological properties, in particular in relation 
with cosmology [11], and some of the drawbacks mentioned above seem to be curable [12]. In the 
perspective of better understanding theories with multiple gravitons, it is interesting to contrast 
theories obtained from a parent theory which is known to be fully consistent, with the parent 
theory itself. This has been done in the past using Kaluza-Klein compactification [13]. Here 
we follow a similar path by comparing simple solutions of theories obtained by discretizing one 
dimension in five dimensional (5D) general relativity to solution of the continuum theory. Namely 
we consider the extra dimension to be given by a one dimensional lattice, to each point of the 
lattice is associated a four dimensional space-time with a metric. We determine the coupling 
of the metrics by the requirement that the continuum limit should be given by five-dimensional 
gravity. The discrete action should respect as much as possible the symmetries of the 5D action, 
otherwise new degrees of freedom appear which are in general ghost-like. The actions we shall 
consider break explicitly the reparametrization invariance along the fifth direction. Although at 
the linear level the theory is free from ghosts, they may appear at the non linear level [4]. 

The paper is organized as follows. In Section 2, we rewrite in a way convenient for discretization 
the 5D Einstein-Hilbcrt action and accordingly recall different ways to parametrize 5D AdS space- 
time as sliced by 4D Minkowski space-time. We then (Section 3) turn to a first naive discretization 
scheme, where the lapse and shift fields (sec below) are omitted. We show in particular that this 
allows one to get a brane world bulk solution with no source brane included. However, this solution 

^and the more recently discussed, and not unrelated, possible appearance of strong coupling [6, 10]. 
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corresponds to a negative tension brane. We show next how one can also recover the positive 
tension configuration, at the price of introducing a fine-tuned brane source in the discretized 
theory. Then we turn to include the lapse in our discretization scheme, and discuss issues related 
to reparametrization invariance along the discretized dimension. In Section 4 we calculate, in the 
model without lapse field, the static gravitational potential between two pointlike particles and 
show that it is 4D at small length scales and becomes 5D at large distances. We then comment 
on the vDVZ discontinuity in this model. In section 5 we discuss cosmological solutions of the 
discretized theory. Our conclusions are collected in Section 6. 



2 Continuum theory: Five dimensional General Relativity 

For the piuposc of discrctizing a space-like dimension, parametrized by coordinate y, it is conve- 
nient to rewrite the 5D Einstein-Hilbert action 



Seh = Mj) J d'X^-giR - 2A), (1) 

using a 4+1 splitting of space time. Above, and in the following, we use expressions with a tilde for 
quantities of the 5D continuum theory, like the 5D metric qab, and upper case Latin letters from 
the beginning of the alphabet. A, B, C... to denote 5D indices, A is the 5D cosmological constant 
and M(5) is the 5D reduced Planck mass. After an integration by parts, (1) is rephrased into 



j d'xdy^H {i? - 2A + K^^K^p (g^^g^^ - g^^g^^) } , (2) 



where iC^i^ is the extrinsic curvature of surfaces Tiy located at constant and we have introduced 
in a standard way^ the lapse J\f, the shift N^, and induced metric g^i, on Tiy whose Ricci scalar 
is denoted by R. Note that here and in the following, lower case Greek letters, ^u, a, ... will be 
denoting 4D indices, i.e. indices parallel to Tiy. The extrinsic curvature is defined by 

K,. = ^ - - D^N,) , (3) 



where is the covariant derivative associated with the induced metric g^j^^ and a prime denotes 
an ordinary derivative with respect to y. The fields jV, A^^ and gij_i, are simply related to the 
components of the 5D metric gAB by 

Qixu = g^iu, (4) 
~gt.y = N^ = g^aN'', (5) 
9yy = Af^ + g^uN''N\ (6) 

The equations of motions derived from the action (2) varying with respect to jV, A^'* and g^" yield 
respectively 

R-2A = K^-RPR^ (7) 



^note however that the surface located at constant y are timehke, unhke in the usual ADM splitting. 
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= D'^K-D'^Ki;: (8) 

+2 [k^^K,, - KK^,) , (9) 

where G^i, is the Einstein tensor for four dimensional metric g^,^, and K is defined hy K = Kf^^g^'^ . 
When vanishes, equations (7) and (9) reduce respectively to 

R-^K ^ ^,9'M9"'9''''-9"'9''% (10) 
G^. + A^M- = ^{D^DM-9,..9"'DpDM) 



+ 



9fj.u 

JP \ 8 



g\9a/39 > + ^9a/39\s9 9 2^ ^ 2j\f j 



+ ^29'ixu{9ai39"'^) + 2_/\/'2^"'^^ 2^/"^^^" 2jV2^^"^^^^"''' 

where G^,^ is the Einstein tensor of the 4D metric (^^j,. We then consider simple solutions of the 
5D Einstein equations, which are slicings of the 5D space-time by 4D Minkowski space-time. We 
write accordingly J\f = J\f{y), gnu{x'^-,y) = ^{y)iliJ.v-i where rj^^ is a 4D Minkoswski metric, and we 
keep Nn — 0. Then equations (8) are identically satisfied and the equations (11) and (10) reduce 
to 

One notes that the above equations are not independent, the first being a consequence of the 
second. This is of course a simple consequence of the y reparametrization invariance of the 5D 
Einstein-Hilbert action (1) and tells us that we cannot solve both for M and Q, out of the equations 
of motion. Indeed by choosing a Gaussian Normal gauge, defined here simply by setting jV' to 1, 
one can get the following solution for the 5D line element 

ds" = -g^dX^'dX^ 



= dy'^ + exp I eyy-^A I dx^dx'^T]^^, (14) 

where we have assumed A to be negative. These solutions are well known parametrization of a 
Poincare patch of AdS^. If one instead puts an absolute value on y, in the above line element. 
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and considers then the metric g]^^ given by 



= V + exp |^e|y|y-^Aj dx'^dx^T?^^, (15) 

one obtains the well known Randall-Sundrum [1] type of metric which are simply given by gluing 
two identical patches of the space-time parametrized by metric g'^^^ along a brane of positive 
(e = —1) or negative (e = +1) tension. This requires a fine tuning between the brane tension a 
and the bulk (negative) cosmological constant A given by 

= -6Mj)A, (16) 

where normalization of the brane tension is set by its action, 

- j d''xdyV^26{y)ea, (17) 

that was added to the Einstein- Hilbert action'' (2) to get (15). Another parametrization 5!^^ of 
the space time described by the line element (14) that we will use in the following, is given by 



;;(3) 

+ Zo{z + zo)v^udx''dx'' , (19) 



ds^ = -gfedX^dX"" (18) 
3 dz'^ 



2A {z + zof 

where zq and Zq are some constants, and z is related to y of equation (14) by Zq[z -\- zq) — 
exp (ey^-2/3A^ 



3 Theories with a discretized extra dimension 

We are now turning to the type of theories which will be of interest in this work. Those can be 
obtained from the Einstein-Hilbert action (2) where one discretizes the continuuous coordinate y 
with a spacing a between two adjacent sites (labelled by an index i). It was shown in reference 
[14] how to obtain such a discretization, maintaining y dependant 4D gauge invariancc on each 
sites. This can be done by the mean of link fields X'^{i,i + l;x) [6], mapping between site i and 
site i + 1, which were explicitly build out of the 5D metric in reference [14]. One can proceed 
as follows, we first note that the y derivatives arising in action (2) only appear in the extrinsic 
curvature Kn^, (as defined in Eq. (3)). The latter can be expressed in term of the Lie derivative 
along the vector field Dy defined by Dy — dy — N^d^, such that 

K,. = ^Cnyg,.. (20) 

''Note that we did not include any Gibbons-Hawking term in the brane action, since this term is cancelled by 
the integration by part to go from (1) to (2) 
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To discretize action (2) in a way which conserve reparametrization invariances on each sites, 
we then simply replace every Lie derivative Coy^ acting on a tensor T{x^,y), by its discrete 
counterpart defined as 

a 

where the action of the transport operator W{i,i + 1) on the components of tensor T is given by 

[W(i, i + m^,]^^,,„,^ (x) = d,,X''^ . . . d,^X^^T,^, {X^{i, i + 1; x)),^,...,^ . (22) 

The transport operators W{i,i + 1) are generated by the shift vector fields [14] and at leading 
order in the transverse lattice spacing a, is given by 

X^'ii, i + l,x)^x'' + aNi'iyi, x) + ^(a^), (23) 

where yi is the y coordinate of site i. If one then makes a gauge choice such that X^^{i, x) — x^, 
one is lead to consider theories of a set of 4D metrics g'^^, and 4D scalar lapse field A/i with actions 
of the form 

S[gi,Xi\ = E^M^^) I d^x^Ni {R{gi) - 2A) - / d'^x^Vigi, g,+,), (24) 

where V{gi, fi'i+i) is an interaction term between the metrics gf^^ and fif^t^, and M(4) is a mass scale 
which sets the coupling scale between the metric on a given site and matter sources that one may 
wish to put on the same site. If one insists in keeping the link with the 5D theory, one should 
verify the equation of motion for X'^ . The latter read (in the gauge = x^) 

- (d.9'^ [9%' - 9'^,) (g^gf - grgf) (25) 

and reduce to equation (8) in the continuum limit. This equation is somehow similar to a Kaluza- 
Klein consistency condition [15, 16]. Note that the index i can be envisioned as labelling theory 
space sites in the spirit of the deconstruction program of Ref . [6] , but theories under investigation 
here can also be considered without an explicit reference to a continuum limit, simply as theories 
of multigravity [3, 5]. In the latter case, one does not have to consider equation (25). We will, for 
most of the cases discussed in this paper, not include any matter fields so that each of the sites i 
will only be considered endowed with a cosmological constant A and will consider various possible 
interaction terms given by 

2 

m 



9^^l) = (^^' - ^M.) (^S' - 9U) {argf - 9^9-') , (26) 

V2{9i,gi+i) = -^V^^{9^;l-9n{9?^-9f){9l.9'a(3-9;a9lp): (27) 
V^{9^, g^^l) = (Ci - an {9at " 9'ap) fc.^f ' ^^5^) , (28) 
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where g^'^ is the inverse metric of g^^^. With such choices of interaction terms, action (24) is a 
simple minded discretization of the 5D pure gravity Einstein- Hilbert action (1), where one has set 
to zero. This can be seen exphcitly using the following identification 

= 1 (29) 

Mj) = Mj^a, (30) 

A - A, (31) 

gU^n = ~9A^^,yi), (32) 
M{xn = M{x^,yi), (33) 

where a is the size of the discretization step along and yi = ia. Note that each of the three poten- 
tial (26), (27) and (28) corresponds to a different way to implement the discretization procedure^ 
outlined in the beginning of this section, depending on whether one discretizes K^i, = Crtygnuf^J^ 
or X'^'^ = -CdX''/'2N. 

3.1 No lapse field and brane space- time with no brane source 

We first set from the beginning the lapse fields Mi to one in the multigravity action S[gi,Mi\ and 
consider solutions to the equations of motion derived from the simplified action for the metrics 
'S'is'tjl]- We wish here to compare these solutions with solutions of the continuum theory 
defined by action (2) and seek solutions of the form 

= ^'V,u, (34) 

with f2* constants. Let us first consider the case of the interaction term Vi. After a straightforward 
calculation, the equation of motion for the metric gi reduces to 

-A = -2 + /, + (/,_!)-', (35) 

with A = 2A/3m^, and fi = fli+i/fli. Equation (35) defines a sequence (see figure 1) which has 
two fixed points (for A obeying A(— 1 -|- A/4) > 0) /+ and /_ given by 




/±= 1-- ±, A --1 . (36) 



These two fixed points are each located on a different side from 1 (for negative A) , which is the 
double root of the fixed point equation for A = 0. Let us first consider the solutions T± defined 
by Vi, /i — f±. This translates into 

= {f±r''' . (37) 



^For simplicity, we have written here only the equation of motion (25) corresponding to the potential Vi 
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We then look at the hmit of this solution when the discretization step a goes to zero. ^From 
equations (29) and (31), one has A = |Aa^, which gives in the limit of small a 




(38) 



leading to, 



2 r 



(]0±exp|±^/--A|/, Ir/^,. (39) 



This matches the metric (14), and the fixed point solutions !F±, in the limit of small a thus each 
corresponds to a discretized version of a Poincare patch of AdS^. 

However equation (35) has more interesting solutions, it can indeed be exactly solved to yield 

n^ = ((/+)^ + /c(/_)^)^, (40) 

where K, is an integration constant, related to /o by 

]C^-l±^. (41) 
/- - /o 

Let us consider the case where /C is positive, this is true if we choose /o to lie in between /_ and 
/+. Note first that, for a fixed value of a different from zero, one can choose any positive value of 
K. One easily finds that the limiting behaviour of the metric on site i, when a is send to zero is 
given by 

^ (iT§'°'^ [f^^^y- ~ ^^^^ 

where T is a constant given by 



It turns out that the solution (42) has exactly the same asymptotics^ as yi goes to ±oo as the 
Randall- Sundrum metric (15) with a negative tension brane placed at y = T. This is a quite 
remarkable feature since no brane has been considered. In particular no fine tuning of the form 

(16) is necessary. 

^Notc however, that the exact expression (40) does not approach the asymptotic form given by (42) uniformly 
with i, i.e., for a fixed value of a, one can make the difference between as given by (40) and the asymptotic form 
read from (42) as large as one wishes going to large value of yi. On a given compact interval though, say of the 
form [—ykiUk] {k being fixed), the convergence is uniform. 
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One might be concerned about knowing to what extent this result depends on the form of the 
discretization scheme. We now investigate this question by searching for solutions of the form (34) 
with V taken of the form 

V^CiVi + C2V2+C3V3, (44) 

by changing we can always assume that Ci + C2 + C3 — 1. This choice also insures that the 
multigravity action S[gi, 1] has the action (2) for a limit when a goes to zero, with the same 
identification as in equations (29-33). It is then straightforward to work out the equations of 
motion for the metrics gi^^ and get for the ansatz (34) the following expression replacing equation 

(35) 

-A = Gi(/0 + G2(/,-i), (45) 



with 



Gi{u) = -4Ci-3C2 + 2 + M(-i + |Ci + ic2 



2 2 2 

u \ 2 2 ' 2 V u 
1/1 3^ 1^\ 1/1 1^ C2 



+r ( -0 + 0^1 + 0^2 ) -7^C2, (46) 



u 



= r(-5 + 5C. + 5C.j + ^:5(5-5C,-5C.J + ^. (47) 



It is easy to see that m = 1 is still a double root to the equation Gi{u) + G2{u) = 0. For A < 0, 
this double root degenerates into two root f±, which are given (for small enough A) by 

/± - 1 = 5± ~ ±a/^. (48) 

Notice that the above expression is independent of Ci and C2, as could have been expected, and 
matches the one found in equation (38). While in the vicinity of these roots the sequence 6i defined 
hy fi — 1 + 5i verifies 

Si^5i^i + Sl-Sl,. (49) 

This leads to the same kind of solution as above: namely a solution with the same asymptotics as 
the Randall-Sundrum solution with a negative tension brane. 

There is an easy way to understand these solutions by comparing the equations of motion 
of the continuum theory (12-13) to the ones of the multigravity theory (35) (or (45)). Indeed, 
in the limit where a goes to zero, equation (35) (or (45)) reduces to equation (12). If one then 
looks for solutions of (12), with jV set to one, one finds that the most general solution is given 



by a linear combination of exp (ky) and exp(—ky) with k given by y — 2A/3. This matches what 
is found in the discretized theory. However, in the continuum theory, the equation of motion 
for J\f, namely equation (13), allows to keep the decreasing or increasing exponential, but not a 
combination of the two. This also enables to understand that we did not find the solution with a 
positive tension brane (which would have been very interesting in many respects). Indeed, a linear 
combination of exponentials exp ±ky (with positive coefficients) is an increasing function for large 
positive y and a decreasing function for large negative y. So that the asymptotic jump of the first 
^/-derivative of the 4D metric across the brane, defined as [^'(+00) —Q'{—oo)]/[Q{+oo) +Q{—oo)], 
is necessarily positive. This is turn means that the brane tension has to be negative, as can be 
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seen from the junction conditions. In the next subsection, we show how to recover the positive 
tension brane solution, introducing a brane source with fine-tuned tension. We then, in the last 
subsection, turn to investigate the discretized theory where the lapse field H is kept, which seems 
to be required if one wants to approach closer solutions to the continuum theory. Let us underline 
here, however, that the solution (40) is a perfectly honest solution of the multigravity theory 
defined by action S'[(?j, 1]. To complete this discussion, let us add that the solution (40), and more 
generally solution of equation (45), can also differ much more dramatically from solutions of the 
5D continuum theory. In particular, the fixed points solutions !F± are very tuned solutions. Indeed 
any departure of /o, say as an "initial" choice from f± leads either to the negative brane-hke 
solution, or to other type of solution where the signature of the 4D metric can change (as can be 
seen from equation (40) with /C < 0). 

3.2 Positive tension brane and fine tuning 

In this subsection (and only here), we include a brane sourcc^^, located at fixed y, in the discretized 
theory considered above. We wish to recover the discretized Randall-Sundrum metric with a pos- 
itive tension brane, and discuss how the fine tuning condition (16) is translated in the discretized 
theory. To do so, we simply add to the action S[gi, 1], considered above with the potential Vi, the 
discretized version of the action (17) which can be simply taken to be 

- j d^'x^/^o^a, (50) 

with a — a, and we have chosen the brane to be located on site i — Q. The only difference with the 
analysis done in the previous subsection is that the cosmological constant A of site i = is now 
replaced by a / M'^^y We now seek to reproduce the Randall-Sundrum metric (15) with e — —l 
in the continuum limit. For i 7^ 0, the equations of motion of g^^^ lead to the same equation as 
(35) that we rewrite 

/^ = i^A(/^-l), (51) 

with Fx{fi-i) = 2 — A — (/j_i)~^. If we consider first the sequence of /j with i > 0, then we need 
this sequence to lead to a decreasing a\ for large i, however the general solution (40) shows this 
is never the case unless we choose fo — f-. Similarly, considering now the sequence of fi, with 
i < —1, we need to choose /_i to be equal to /+ in order to get an increasing function of y, as y 
becomes large and negative. In order to reproduce the Randall-Sundrum solution, with a positive 
tension brane, one thus needs to jump from /_i = /_|_ to /o = /_. However f_i and /o are related 
by the g'^^^ equation of motion leading to 

fo = -^A+2o-/3m2M2^j (/-l) ■ (52) 

^°The term "initial" might be somehow misleading here, since the sequence fi is also continued in the "past", 
i ^ — oo, from the "initial" value /q. 

^^Namely, in contrast to the previous cases, we are now explicitely breaking the i translation invariance by 
allowing the " cosmological constant" (or tension) of a given site (or brane) to differ from the one of its neighbours. 
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Demanding then /_ = F^+2a/3m'^M^^^{f+) leads to a fine tuned value for a given by (see figure 1) 



(53) 



This in turn, leads back in the small a limit to the constraint (16). 




3.3 Lapse field and reparametrization invariance 

We then turn to another possible discretization of the action (2) keeping the lapse field on all 
sites. Starting with V — Vi, the equations of motion now read 

(54) 

m'^'-^-M'-t:)' '''' 

where = ±1. Fixing independently of i, and replacing by its expression as a function of 
fi given by equation (54) into equation (55), one finds that fi verifies equation (35) with A = 0, 
namely 

-2 + /i + (/i-i)-' = 0. (56) 
This can be straightforwardly solved to yield 

9^ = Io{i + io)Vi^u, (57) 

^fi = J^.y (58) 
-X{i + lo) 

with io and Jq some constants. This exactly reproduces the metric g)^^, of equation (18) sampled 
at successive sites i parametrized by their coordinate Zi = ai along the z direction. Indeed (jf^^^ and 
Mi as defined by equations (57) and (58) are given by gf^^ = g'^^^y^Zi^ x^) and Mf = g^^J{zi, x^) with 
the identifications Zq — IqQ and Zq — lo/a. This means in particular that our discretization scheme 
reproduces, when the discretization step is send to zero, solutions of the 5D continuum theory 
in a given gauge. This is of course not a surprise, since our discretized action S[gi^Mi\ explicitly 
breaks reparametrization invariance along the y direction. It can also be seen by noticing that, 
in contrast to the continuum case, the equation of motion of the metric (55) cannot be obtained 
from the equation of motion for the lapse field (54). This very fact was indeed used to solve for 
Mi and (yf^^. However, note that both equations reduce to equations (12) and (13) in the leading 
order in a. So that once the solution for Mi (respectively g''^^,) of equations (54-55) is known, which 
can be considered from the point of view of the continuum theory as setting the gauge, as far as 
y reparametrization is concerned, then the solution for the gf^j, (respectively Mi) agrees with the 
solution of the continuum in the small a limit. 

Another property of the solution (57,58) is that the parameter allows one to avoid the 
potential singularity which exists in the continuum. In fact if Iq is not an integer than the scale 
factor Oj does not vanish for all i. This is a generic feature of discretized models. 
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Figure 1: Depiction of the resolution of equation (35) without (subsection (3.1)) and with (subsec- 
tion (3.2)) a brane source included. Equation (35) defines a sequence which for A < and small, is 
characterized by the upper solid curve. The latter curve is representing the function F\{u) which 
defines the sequence by the recursion relation fi = Fx(/i_i). The solid oblique line is defined by 
the function u — v. One sees that, when no brane is included, one can only flow (black arrow) from 
the fixed point /_ to the fixed point /+, for an "initial" value /o which lies in between the fixed 
points. This corresponds to the solution (40) and to the negative tension brane discrctized bulk. 
When a brane source is included, one can go "backward" (white arrows and double dashed lines) 
from the fixed point /+ to the fixed point /_ at a given value of the index i corresponding to the 
brane. For this to happen, the brane tension has to be tuned in order for the curve representing 
-^A+2(T/3m2M2^j) the dashed bottom curve, to go through the point {u — f+,v — /_). The solution 
is then the discretized Randall-Sundrum metric with a positive tension brane. The dotted curve 
represents the function Fq. It is tangent to the u = v line at the point (1,1). It also corresponds 
to the sequence which solution is obtained in equation (57) of subsection 3.3. 



One can ask whether one can restore the arbitrariness in Mi by changing the discretization 
scheme. It turns out indeed to be possible choosing the interaction potential V of S[gi,N'i\ to 
depend on metrics on three adjacent sites, and reading 

V{gi-,, 9.. Qi^,) = -^v^ (9i^ - 97^ - ^r/) {9T9f - ^r^f ) ■ (59) 

With this interaction term, the action S\gi^Hi\ also agrees with action (2) in the continuum limit. 
The equations for and Mi now read 

If one chooses the same value of for all i, then the first equation is a consequence of the second. 
It is thus possible to "choose a gauge" by fixing Mi and then to find Jlj by solving (61), the 
resulting solution will automatically satisfy (60). Let us choose, for instance. Mi — \ — ei^ then 
(61) can be readily solved to give a solution depending on two parameters 

= q;xY + /3xL, (62) 

where 

x± = ± Vl- A. (63) 

Notice that a;_ is negative so that the solution with a = has Jl' with alternating signs: neighbor- 
ing spacetimes i and i + 1 have opposite signatures. This behavior has no continuum counterpart 
where a singularity is reached before the signature change takes place. When /9 = the solution 
coincides, to the first order in a with the solution we found prcviouscly in (37) and so it tends to 
the AdS^ metric in the continuum. Notice that had we started to solve (60) with Mi = \ than 
we would have obtained the analog of the more general solution (40) but equation (61) would not 
have been satisfied. 

To summarize, the inclusion of Mi allows to recover just the continuum solutions if we do not 
insist on having an Mi undetermined by the equations of motion^^; and if we do insist, we get 
extra solutions with no continuum counterpart. 



— eiTTi 



a 



+ 



VtiMi-i 



(60) 
(61) 



4 Effective 4D gravity 

We turn now to determine the gravitational potential without the lapse field, and consider the 
discretisation scheme given in the previous subsection. We assume here that the cosmological 
constant vanishes and so we can perturb around fiat space-time. To quadratic order in the metric 
perturbation /i^^, the potential (59) reads 

2 

V{K.,, h„ h,+,) = - h'-j) - {V'^'V''^ - V'^V''^) ■ (64) 

^^and we do only consider solutions which are continuous in the continuum limit 
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In order to diagonalize this interaction, we define h^i, by 

V(^",^)=E^M^(^")e'"'- (65) 

n 

The quadratic action is then an integral over 6 of Pauli-Fierz actions with a continuous mass 
spectrum given by 



m%9) = — —. (66) 



a" 

And the gravitational potential ^j(r), between two unit masses separated by r and placed at 
sites i and i + j, can be obtained summing over Pauli-Fierz propagators. The outcome of the 
calculation can also be simply understood from discrctizing the Laplacian equation with the same 
discretization scheme as above. The discretized equation reads 

a=*,+ (5i±i+^^ii^ = 4.G.^(r)^,„, (67) 

where we have reintroduced the Newton constant Gn- The Fourier transform ^, as defined above, 
verifies 

d^^{r, 9) - m^{e)^{r, 9) = 47rG'^5(r). (68) 

Notice that the mass spectrum is bounded from above by the inverse lattice spacing a^^. A 
continuous mass spectrum is reminiscent of the infinite dimensional models of [19, 20]. The 
gravitational potential can now be readily put in the form 

^.(^) ^ _^ £ ^e-S-^^osj^. (69) 

When r << a then the integral can be approximated by Sj^ and the potential reduces to the 4D 
Newtonian potential 

^Ar) = -^{ko + oC-)), r«a. (70) 
When r » a, then the integral can be approached by 

r^'^e-^^osj^+ r -e-S(-^)cosj^ = -^-^(l + (-l)^ + 0(e^)), (71) 

Jo TT Jn/2 TT TTO + O^J^ 

SO that the gravitational potential is of the form of a 5D potential 

= -^^7^2^ (72) 
if j is even, with the 5D gravitational constant being given by 

G, = (73) 
na 

When j is odd the gravitational potential is exponentially small. Gravity is thus four dimensional 
at small length scales and five dimensional at large scales. Choosing a very large (of the order of 
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the Hubble scale) allows a very simple modification of gravity at large scales in the spirit of the 
models of references [18, 19, 20, 21]. Note that we could have started from a finite number N of 
sites. This corresponds to an extra dimension which is compact with a length scale it! = Na. The 
graviton spectrum would have been discrete with a spacing of order 1/A^ and would have remained 
bounded. In this case, the four dimensional regime is obtained both for small scales r << a and 
large scales r » R, the intermediate scales being five dimensional a « r « R. The qualitative 
features of this potential are not sensitive to the particular form of the discretization scheme we 
have used. On the other hand, it is well known that the tensorial structure of the propagator 
of massive spin 2 fields differs dramatically from the massless one. This leads to the vDVZ 
discontinuity at the linearized level which is manifested by, e.g., an order one difference in light 
bending. In this respect, when R is infinite we expect the discontinuity to be present since we 
have a continuous spectrum of massive gravitons. This is similar to the brane models with an 
extra infinite dimension [18, 19, 20, 21]. When R is finite, however, the spectrum is discrete and 
there is no discontinuity. 



5 Cosmological solutions 



The discussion of preceeding subsections enables to find easily some cosmological solutions for the 
type of models considered here. It is the purpose of this last section to discuss those solutions. 
We look for solutions where the line element of site i takes the form 

dsl = -Bl{t)de + A]{t)dx'. (74) 

The equations of motion derived from action (24) (where we consider Hi as a non dynamical field 
set to 1) with the potential Vi read 

A? ' 2a2 r ^ 25? 2AW} 2A| ^ 2^ J 

+ M(-2i?2p., (75) 




A "Ai ' "AiBi ' 2a^Af \ 2Bf 2Af 2 ' 2Bf 



where we have included the possibility for each site i to host matter with energy density pi and 
pressure pi (that is to say, matter, in the form of a perfect fiuid, is assumed to be confined on each 
site) . We then look for factorized solutions of the form 

A(t) = a,,A'{t), 

B^it) = (3,B'{t), (77) 

where and f3i arc timc-idependent, we also set B to one by a suitable time redefinition. We 
then further simplify the problem by letting ctj = and look for solution where the metric on 
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each site is de Sitter, and the sites are only endowed with a cosmo logical constant a. This means 
Pi obeyes pi = —pi = a and does not depend on i. In this case, the equations of motion reduce to 



— (a^+i + 2q;0 = 3/^2 _ (a + M(^)V) . (78) 

These equations can easily be solved. Let us first discuss the case where (A + M^a) does not 

vanish. In this case, defining 7 by 7 = 3H'^/ (Ji + M^^^a^, the sequence = (oj+i — 'y)/{o(i — 7) 

obeyes equation (35) with A given by A = 2(A + M^^a)a'^ /3. So that the solutions are readily 

obtained from those of section 3. In the case where (A + M^^a) vanishes, on the other hand, the 
solution to equation (78) is given by 

ai — ao + agi + H'^aH^ , (79) 

where and do are some arbitrary constants. Interestingly, this solution has de Sitter 4D space- 
time on each sites, but vanishing cosmological constants. This is in analogy with what is happening 
in the brane-induced gravity model [22, 11] as is the behavior of the Newtonian gravitationnal 
potential that was discussed in the preceeding section. However in contrast to brane-induced 
gravity, the curvature of de Sitter space-time is here arbitrary. Note that, had we kept the Xf 
fields as a dynamical variables, one should have also considered their equations of motion. The 
latter, for the ansatz (77), reduce to 



11 2^ _ 3 p±i _ ^:!±i = 0. (80) 

This equation is satisfied for our choice ai = (3i. However, in a similar way as was discussed 
previously, solution (79) does not verify (for H ^ 0) the equations of motion for A/i which reads 
(forM = l) 

^/o^V (81) 

Indeed in the continuum theory, there is no solution equivalent to (79) with arbitrary ao and (5o; 
rather one finds a related solution describing a shcing of 5D Minkowski space-time by 4D de Sitter 
space-times which reads 

ds^ = (1 + Hyf (-dt" + e'^'df) + dy\ (82) 

This solution is conveniently used to described the space-time on one side of a domain-wall, which 
entire space time is given by the above solution where y is replaced by |y| [23]. 
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6 Conclusions 



In this paper, we have investigated some exact solutions of 5D General Relativity with a discrete 
fifth dimension, and compared them to solutions of the continuum theory. Depending on the 
discretization scheme used, we have shown that some of the solutions of the discrete theory exactly 
match those of the continuum, while others do not. In general, the discretization explicitly breaks 
reparametrization along the discrete dimension. This gets reflected in the fact that for solution we 
have considered, corresponding to slicing of AdS^ by Minkowski space, the equations of motion 
of the discrete theory enable in general to determine both the lapse field and the 4D metric. 
This is to be contrasted with the continuum theory, where the lapse field cannot be determined 
by the equations of motion. However, we have shown that one can find a given discretization 
scheme in which it remains undetermined. More importantly, some of the solutions of the discrete 
theory exibit very dramatic differences with those of the continuum. As shown here, one can 
find signature change of the 4D metric but also avoidance of singularities which would be present 
in the continuum. In addition we have also found a brane world looking bulk space time, with 
no brane source in the equation of motion and no fine tuning. This solution would correspond 
in the continuum theory to a negative tension brane. We investigated the static gravitational 
potential and found it realistic when the lattice spacing is very large. Last, we also discussed some 
cosmological solutions. 

We can look at our results from two different perspectives. On one hand, they exemplify the 
difficulties which arise upon deconstruction of gravity even at the classical level: the neglect of the 
lapse fields leads to spurious solutions, while their inclusion only partially solves this problem, as 
we showed in section 3. On the other hand, they point out interesting directions in multigravity 
theories: they allow, as we showed, simple modifications of gravity at large scales and brane-like 
solutions with no branes. One can speculate on the possibility to reproduce these results, in 
some more complex multigravity, without the various drawbacks mentioned in the introduction. 
It would also be interesting to find a solution corresponding to the discretization of RS type - 
gravity localizing - metric [1], and at the same time avoiding fine tuning conditions. 
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